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Several heavy-fermion metals display a quantum phase transition from an antiferromagnetic metal 
to a heavy Fermi liquid. In some materials, however, recent experiments seem to find that the 
heavy Fermi liquid phase can be directly tuned into a non-Fermi liquid phase without apparent 
magnetic order. We analyze a candidate state for this scenario where the local moment system 
forms a spin liquid with gapless fermionic excitations. We discuss the thermal conductivity and spin 
susceptibility of this fractionalized state both in two and, in particular, three spatial dimensions for 
different temperature regimes. We derive a variational functional for the thermal conductivity and 
solve it with a variational ansatz dictated by Keldysh formalism. In sufficiently clean samples and 
for an appropriate temperature window, we find that thermal transport is dominated by the spinon 
contribution which can be detected by a characteristic maximum in the Wiedemann-Franz ratio. 
For the spin susceptibility, the conduction electron Pauli paramagnetism is much smaller than the 
spinon contribution whose temperature dependence in three dimensions is logarithmically enhanced 
as compared to the Fermi liquid result. 



I. INTRODUCTION 

In many heavy fermion systems, the heavy Fermi liq- 
uid phase seems to disappear at a zero temperature 
phase transition upon tuning a certain external control 
parameter.^ Since the discovery of such transitions, it has 
remained an open problem what precisely happens at the 
underlying quantum critical point (QCP). One peculiar- 
ity is that often, a simultaneous occurence of magnetic 
order accompanies the breakdown of the heavy Fermi 
liquid at the QCP. An important question in this con- 
text has been whether the nature of magnetic order is 
local or itinerant close to the QCP where different ma- 
terials appear to provide different indication.'^' Of partic- 
ular interest are materials where the heavy quasiparti- 
cles seem to lose their integrity across the QCP.^^l In this 
case, it is challenging to understand how magnetism can 
emerge out of the heavy Fermi liquid, since no conven- 
tional mechanism such as a spin-density wave instability 
is applicable.^ 

In order to concentrate on the role of Kondo screening 
across the transition, theoretical concepts have been dis- 
cussed that ignore long-range magnetic fluctuations at 
the quantum phase transition, such that the local mo- 
ments enter a spin liquid state upon tuning the system 
across the QCP. Although this fractionalized Fermi liq- 
uid has been shown to be a stable phase of matter,*^ 
in heavy Fermion materials, it was long believed that 
ordered ground states are more natural to occur. Re- 
cent measurements on chemically substituted YbRh2Si2 
suggest a different situation.^' Upon partially replacing 
Rh by Ir, magnetic order and the heavy Fermi liquid 
state appear to separate and a non-Fermi liquid state 
without any apparent magnetic order extends over a fi- 
nite range of magnetic field strength. Similar observa- 
tions of non-Fermi liquid behavior without indication 



for a magnetic QCP have been reported in other Yb- 
based three-dimensional heavy fermion systems, includ- 
ing Ge-substituted YbRhaSiJ^and /3-YbAlB4.'^ Another 
relevant material is the geometrically frustrated Kondo 
lattice PrsIrsOy.* 

In this paper, we consider a fractionalized Fermi liq- 
uid state as proposed by Senthil et al. as a possible 
candidate state for such phases.'sEl There, the Kondo- 
Heisenberg lattice model is used as its microscopic start- 
ing point. In a large-N treatment, such a model naturally 
describes the breakdown of Kondo screening at a QCP 
where the heavy Fermi liquid is destroyed. This break- 
down is described by the hybridization mean field be- 
tween local moments and conduction electrons. Within 
a large-N theory, the hybridization amplitude vanishes 
at the QCP, as well as the Kondo temperature derived 
from the hybridization mean field.l^Upon including gauge 
field fluctuations into this model, the transition becomes 
a deconfining transition of a U(l) gauge field such that 
the local moment system forms deconfined spinon exci- 
tations. The resultant state has been coined fractional- 
ized Fermi liquid, abbreviated by FL* . The theory shows 
several appealing features that are consistent with exper- 
imental observations.!^ One particular result is a loga- 
rithmic divergence of the specific heat coefficient at the 
QCP, as measured in CeCue-ajAua;.!^ Another prediction 
is a jump of both electrical conductivity and Hall coeffi- 
cient across the QCP. "'^^ The latter has been inferred from 
measurements on YbRh2Si2 compounds ,1^^ while differ- 
ing interpretations of the findings in these materials are 
also possible. ^'^ Whereas several subsequent works have 
analyzed properties of thi s model at the quantum crit- 
ical point in great detail ,1^^3111111 not much attention 
has been assigned to explicit predictions of thermody- 
namic properties of the FL* phase that are measurable 
in experiment. Although recent experimental efforts have 
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promised novel possibilities to explore such phases inma- 
terials based on the parent compound YbRh2Si2P^ it 
still remains difficult to identify such phases in terms 
of concrete measurements: in addition to the fact that 
spin liquid states do not break any lattice symmetries, 
fermionic particles in a U(l) spin liquid state also obey 
the kinematic constraint of not coupling to the external 
electrical field. The existence of such degrees of freedom 
hence cannot be inferred from their contribution to the 
electrical conductivity or the Hall effect. 

The determination of measurable signatures of a frac- 
tionalized Fermi liquid is the main task of our paper. 
The angle from which we intend to tackle this problem is 
that despite of not carrying electrical charge, fermionic 
spinous are still able to carry heat and to respond to 
external magnetic fields. These properties can be used 
to analyze signatures of fractionalized phases and help 
to either substantiate its existence or non-existence in 
forthcoming experiments. In this paper, we discuss the 
thermal conductivity and the spin susceptibility for the 
fractionalized Fermi liquid candidate state proposed in 
Ref. 4j Our results are based on a low energy effective 
theory for a U(l) fermionic spin liquid state. As shown 
previously,'^the fermionic spinon excitations do not count 
to the Fermi volume and can be treated independently 
from the conduction electron degrees of freedom. The 
stability of such a state beyond the mean filed limit is 
still under debate in two spatial dimensions, while its ex- 
istence is established in d > 3 which is the main focus 
of the paper: building up on previous works for the two- 
dimensional case, we particularly investigate the three- 
dimensional setup and compare with previous results for 
the 2d case. 

Starting from the low-energy effective theory, we derive 
the Keldysh equations of motion for the non-equilibrium 
spinon Green's functions in presence of a thermal gra- 
dient. We linearize these equations of motion and argue 
that they can be rewritten in form of a variational princi- 
ple that allows to derive the thermal conductivity from a 
suitable variational ansatz for the non-equilibrium spinon 
propagator. Collision processes in our description origi- 
nate both from impurities and collisions with U(l) gauge 
bosons. In a temperature regime above the impurity- 
dominated regime of the thermal conductivity, we argue 
that the spinous give rise to the dominant contribution 
to the heat current as compared to the heat current car- 
ried by conduction electrons, phonons, or gauge bosons. 
We will explain that the spinon contribution leads to a 
characteristic maximum in the temperature dependence 
of the Wiedemann-Franz ratio at a temperature that is 
strongly dependent on the impurity concentration in the 
sample. 

Our discussion of the spin susceptibility is based on 
the free energy contribution of the combined system of 
spinous and gauge bosons. Based on this quantity, we 
obtain that the temperature dependence of the spinous 
is logarithmically enhanced as compared to the temper- 
ature dependence of conduction electrons in a Fermi liq- 



uid scenario. We find the zero temperature paramag- 
netic susceptibility contribution of the conduction elec- 
trons and spinous to be larger than the conduction elec- 
tron Pauli susceptibility and thus the low-temperature 
susceptibility promises to be a suitable quantity to iden- 
tify fermionic spinon excitations. 

The theory of a fractionalized Fermi liquid has been 
originally developed in the conte xt o f unconventional 
QCPs in heavy fermion materials.'^ Lateron, several 
more detailed discussions of the quantum critical regime 
associated with this transition appeared. I^ESIIll par- 
ticular, predictions have been made for a diverging 
Griineisen ratio and a violation of the Wiedemann- 
Franz ratio.l^^E^ Furthermore, spatially inhomogeneous 
solutions in the heavy Fermi liquid phase have been 
discussed,'^^ and it has been predicted that volume col- 
lapse transitions are likely to occur near such quantum 
phase transitions in sufficiently soft materials.C^ Neither 
such first order transitions nor inhomogeneous solutions 
in the Kondo screened phase could be uniquely related to 
the scenario within the FL* phase. Instead, we attempt 
to approach the problem directly from the FL* trial state 
and derive its signatures for some directly measurable ob- 
servables. 

It is frequently stated that the fermionic U(l) spin liq- 
uid has a stable low-energy fixed point in the limit of large 
number of fermion flavors. In two dimensions, the exis- 
tence of such a limit has been questioned recently,'^^! due 
to the existence of infinitely many corrections to the low 
energy fermion propagator from planar diagrams. The 
subsequent proposal of an expansion in the Fermi surface 
genus has been challenged by Metlitski and Sachdev.l^ 
In our work we will not elaborate on the existence of a 
large-N limit and just take the previous 2d results as valid 
and given, assuming that the corrections to the one loop 
results in 2d are small for finite N. In three spatial di- 
mensions there exist no singular corrections from higher 
loop contributions which is why we can disregard these 
problems in our calculations. 

Since self-energy corrections due to the fluctuating 
gauge field are singular, transport properties cannot be 
described by assuming well-defined quasiparticles. Kim 
et al. circumvented this problem by formulating a quan- 
tum Boltzmann equation (QBE) for a generalized distri- 
bution function.'^i' This approach has been implemented 
in form of a variational principle by Lee and Nave in 
order to discuss thermal transport properties in two 
dimensions. "- Properties of the magnetic susceptibility of 
two-dimensional fermions coupled to a gauge field have 
been discussed by Lee and Nave.— Similarly, loffe and 
Kalmeyer have discussed these issues for bosons coupled 
to a gauge field.'^^l Building up on these previous works, 
we will in particular compare the contributions from the 
spinous and conduction electrons to thermal conductiv- 
ity and magnetic susceptibility, and mainly discuss three- 
dimensional candidate materials. 

The remainder of the paper is organized as follows. In 
Section [nj we introduce some technical details of the low 
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energy effective theory we will use to describe the FL* 
phase. In Section we discuss several basic concepts 
that are required for a transport theory involving the 
effective low energy description. In Section |IV[ we use 
these aspects to formulate a quantum Boltzmann equa- 
tion approach for systems in presence of thermal gradi- 
ents and derive a solution to this transport equation by 
using a variational principle. Furthermore, we discuss 
the different contributions to the total thermal conduc- 
tivity in different temperature regimes. In Section |Vj 
we discuss properties of the spin susceptibility stemming 
from its spinon contribution. We finally conclude in Sec- 



tion VI that heat transport and magnetic susceptibility 
are promising measurements to reveal thermodynamic 
signatures of the fractional Fermi liquid. 



II. MODEL 

A microscopic starting point for an FL* phase is a 
Kondo-Heisenberg model, both in two and three spatial 
dimensions. Besides the limits of one or infinite spatial 
dimensions,!^ a detailed solution of this model is still out 
of reach, and the existence of a fractionalized phase in 
this model has been only established in the limit of large- 
N. Very recently, an application of the ADS-CFT corre- 
spondence has been used to describe such FL* phases 
in a novel way, starting from the more general Ander- 
son lattice niodel.'^ In this work, we will not add to the 
discussion of microscopic origins of the FL* phases, but 
concentrate rather on its thermodynamic properties. We 
start with the fermionic local moment representation 
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where F" are the generators of the SU{N) group in 
the fundamental representation. The fermions /^^ obey 
canonical anticommutation relations as well as the local 
constraint 
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This constraint is enforced by a U(l) gauge field which 
can be minimally coupled to the spinous. In all following 
calculations, we will set iV = 2, with only slight mod- 
ifications required for general values of iV. In general, 
the fermionic degree of freedom can be complicated by 
its coupling to the conduction electrons, described by 
the bosonic hybridization field fe| — 'Y^^ fia'^ia- We will 
call the particle excitations of this field "hybridization 
bosons" in the following. At zero temperature, the frac- 
tionalized Fermi liquid state is reached technically by 
tuning the chemical potential of the hybridization bo- 
son bi according to fib ^ Jk — Jkc, see Fig. [T] The 
Kondo coupling Jj^ itself can be tuned by variables as 
pressure or chemical composition. The size of the mass 
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FIG. 1: Crossover phase diagram for the Kondo breakdown 
transition. In this paper, we analyze transport properties in 
the (7(1) FL* phase. Upon increasing temperature in tliis 
phase, transport properties are expected to show a crossover 
due to the fluctuating bosonic hybridization field bl in the 
quantum critical region. Figure taken from Ref. |4l 



gap fib is related to a crossover temperature scale T* 
for the FL* phase. In three dimensions, T* follows 
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Below temperatures of 



order T*, the number of hybridization bosons is exponen- 
tially small as a function of temperature, and the system 
of spinous is then described by a U(l) gauge theory which 
is decoupled from the conduction electrons. In this case 
to which we will constrain our attention, the effective 
action assumes the form 
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D{q,uJny^{Sij - ^^)ai(q, w„)aj(-q, -a;„) . 

(3) 

The conduction electron degrees freedom Cko- are en- 
tirely described by their excitation energies ek- The 
fermionic spinous /ko- are coupled to a U(l) gauge field 
with transverse modes a.i and longitudinal mode a,- . The 
gauge field dynamics is entirely generated by the spinon 
matter field, with the retarded propagator of the trans- 
verse fiuctuations at long wavelengths given by D{q, w) = 
{iT{uj/q)+Xdq^)^^- Here, the dynamics of the longitudi- 
nal gauge field mode a,- is neglected, since its fiuctuations 
are short-ranged. Assuming that the spinous have mass 
nif and Fermi wavevector kp, the parameters F = ^ 

and Yd — yield in 3d, while we have Yd — — 

and F = A:f in 2d. The action Sf+Sa has been discussed 
in many different contexts. 1^ Here, we consider it as a 
starting point in order to analyze physical observables of 
the FL* phase. 
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III. GENERAL TRANSPORT PROPERTIES 

A. Linear response 

We are interested in the response of the total electrical 
current density J and the total heat current density Q 
to an applied electrical field E and a temperature gradi- 
ent VT, weak enough in order to linearize J and Q in 
the applied fields. These properties may be formulated 
individually for each type of particle, which we label by 
a = b,Cj f for bosons, conduction electrons and spinous, 
respectively. The transport coefficients for these individ- 
ual types of particles are defined by the relations 



-VT 

T 

-VT 
T 
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where the Onsager relation for thermoelectrical trans- 
port coefficients has been used.^'' The field coupling to 
the electrical transport coefficient Koa can originate both 
from an electrical field or a gradient in the chemical po- 
tential. Precise definitions for the current densities J 
and Q will appear in concrete calculations in Section [IV] 
The thermal conductivity k is defined by the relation 
Q = — kVT with the boundary condition J = 0. This 
implies that in general, the total thermal conductivity 
K will depend on all of the coefficients entering Eq. Q . 
This complicated dependence can be slightly simplified 
by rewriting the thermal conductivity in form of a varia- 
tional functional as described by Ziman,!^ such that not 
all of the 9 coefficients K{oA.2}a need to be evaluated 
explicitly. We will introduce this approach in explicit 
calculations later on. Moreover, any concrete calcula- 
tion becomes simplified considerably if the contribution 
of conduction electron and spinon degrees of freedom to 
the total thermal conductivity can be calculated sepa- 
rately. This case is not justified in general, since Eq. Q 
couples different electrical fields Eq to each type of par- 
ticle. More precisely, corrections to the approximation 
of independent spinon and conduction electron thermal 
currents arise for temperatures above or the same size as 
the crossover temperature T* related to the boson mass 
gap fib- A formal discussion of transport properties in 
this general situation can be given in terms of t he co m- 
position rules first discussed by loffe and Larkin.E^^ 



loffe-Larkin rules 



The fact that the boson variable 6| is given by 6| — 

J2cr fia'^ia immediately implies the relation /if, = M/ ~ 
between the chemical potentials of boson and fermions, 
while the electrical fields felt by the particles are 



Here, the electrical field e is the fictitious field strength 
felt by the gauge charge of a spinon in presence of a finite 
external electrical field E. Since both bosons and spinons 
couple to the internal gauge field a, the gauge current 
densities will fulfill the constraint Jb + J/ = 0. The 
physical electrical current density is given by J = Jc+Jfc. 
From these constraints, one can readily derive the compo- 
sition rules stated in Ref. 1161 representing extensions to 
the composition rules of loffe and Larkin.Eil For the total 
thermal conductivity kj, these rules imply the fornPsl 
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Here, Kot is the total electrical conductivity and ku the 
total thermoelectric conductivity. The correction terms 

^ ~ become exponentially small as a function 

of temperature below the crossover temperature scale 
following T* ^ iMbP^'^ in three dimensions,!^ since the 
number of bosons decreases exponentially as function of 
temperature in this regime. In the same way, Kot and 
Kit reduce to the conduction electron contributions kqc 
and Kic, since the spinons carry no physical electrical 
charge. It might be possible instead to measure the spin 
conductivity of the spinons, which shows a temperature 
dependence identical to the gauge charge conductivity. 
Its temperature-dependence behaves as ~ y5/3 three 
spatial dimensions.^ We do not elaborate on the spin 
conductivity as it is a rather complicated quantity to be 
measured experimentally. 

Above the crossover scale T*, the presence of hy- 
bridization bosons provides both a strong source of scat- 
tering and a sizable channel for heat conduction, and 
we cannot neglect their contribution to transport observ- 
ables. Throughout the following calculations, we assume 
that T < T* and neglect therefore the presence of the 
hybridization boson degrees of freedom. Therefore, we 
can assume that the conduction electron system shows 
conventional Fermi liquid properties, while the spinon 
system will be responsible for non-Fermi liquid proper- 
ties. 



C. Spinon self-energy 

In order to analyze scattering processes of spinons on 
gauge field fluctuations, it is convenient to consider the 
transverse gauge, V • a = 0. In this case, the scalar 
(ao) and transverse parts of the gauge field fluctuations 
decouple, oq corresponds to a density-density response 
function and does not lead to singular corrections to the 
spinon dynamics.l^Sl It is therefore possible to concen- 
trate on the transverse fluctuations. The spinon self- 
energy due two single gauge boson exchan ge processes 
has been discussed by many authors .l^^l^ Close to the 
spinon Fermi surface, the spinon self-energy corrections 
are 
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in 2d and 



I]a(k, uj) — Xiuj In ( — I + 1X210 



(8) 



in 3d.l^ Here, Ai and A2 are constants that depend on 
details of the spinon Fermi surface. Due to these self- 
energy corrections, the quasiparticle lifetime is not much 
longer than its inverse energy and several non-Fermi liq- 
uid properties result both in two and in three spatial 
dimensions!^ An additional complication arises at finite 
temperatures, since Ea(k,a;) is divergent in d < 3 for 
UJ < T.^^ In order to describe transport properties of the 
coupled system of spinous and gauge field, Kim et al. 
formulated a quantum Boltzmann equation approach in 
which the low energy fluctuations of the gauge field arc 
split off from the collision integral in order to avoid the 
divergent self-energyPH 

To cure the finite temperature divergence of the single- 
particle self-energy, vertex corrections are of central im- 
portance. The divergence mentioned above is caused by 
the fact that the single-particle spinon Green's function 
is not invariant under a U(l) gauge transformation.'^ 
In order to calculate a gauge invariant quantity such as 
thermal conductivity, it is also possible to include the 
low-energy gauge fiuctuations into a quantum Boltzmann 
equation, since the divergences are expected to cancel in 
the gauge invariant conductivities. This situation occurs 
in our following calculations. 

At low frequencies the spinous are also influenced by 
impurity scattering for which we assume potential scat- 
tering with potential strength V{r). Any type of im- 
purity with short-ranged Coulomb potential will lead to 
such scattering processes. We assume in the following 
that the scattering potential is independent of tempera- 
ture, as well as a temperature independent concentration 
and a Gaussian distribution of impurities in the bulk of 
the sample. Modifications to the impurity potential are 
expected to occur only at temperatures of order the ion- 
ization energy required to change the impurity charge, 
which is considerably higher than room temperature. To 
lowest order in the impurity concentration n^, such im- 
purity scattering processes are described h^p^ 



Simp(k, Cj) = UiTkk 

with the T-matrix given by the solution to 
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Here, ~ p^/(2m/) — /i is the spinon dispersion and 
V{p) is the Fourier transform of V{r). In the follow- 
ing, we neglect interference processes between scatter- 
ing on gauge bosons and impurities, such that the total 
self-energy is given by E(k,cj) = Einip(k,a;) -I- Eo(k,a;). 
This approximation is valid if the mean free path de- 
rived from one of these two collision processes is much 
larger than that of the second collision process. Since 



the mean free path is given by the spinon Fermi ve- 
locity times the spinon lifetime VF/'E{'kF,uj), it is suffi- 
cient to compare the size of the two self-energy contribu- 
tions Simp(k,w) and Ea(k,w) for this purpose. Setting 
them equal defines a characteristic temperature which 
we denote by Ti. First, the neglect of interference 
processes is therefore valid in the temperature regime 
Ti <^T <T* where scattering on gauge field fluctuations 
dominates, such that E(k, w) ~ ^^(kjti;). The temper- 
ature scale T* is the crossover temperature for the FL* 
phase boundary we introduced above. Second, our ap- 
proximation is also valid for temperatures T <^Ti where 
impurity scattering dominates the self energy such that 
E(k,w) ~ niTkpkF('^ = 0)- The temperature scale Ti is 
set by equating the imaginary parts EJ^p(kF,aj = 0) = 
E"(kF,w) and considering fiw ~ fc^T. We abbreviate 
the impurity self-energy as the spinon scattering rate tJ^ 

such that Einip(kF, lo = = tJ^ . For the constant A2 in 

Eqs ([7]) and ([8|, wc use the estimate A2 ~ {tpY'^ in 2d.l^ 
In 3d, an estimate for metals is given by A2 — 10~^,l^ al- 
though the gauge field coupling of the spinous can differ 
to some extend from the case of electrons in a metal. In 
this way, we obtain 



3/2 



(11) 



in 2d and 



ksTx ~ 10^ 



(12) 



in 3d. We thus conclude that the temperature window 
where the neglect of interference processes is not justi- 
fied is proportional to t^^ ^'^^ ~ nf '^^^'^ and therefore 
shrinks rapidly with decreasing impurity concentration. 
Concise numbers for the temperature Ti can obtained 
by estimating Tf from the residual conductivity expres- 
sion (35), which itself could be estimated from the con- 



duction electron residual conductivity. Another way to 
check the justification of the neglect of interference pro- 
cesses is to compare the mean spacing between impurities 
and the scattering length on gauge bosons. The latter 
will shrink as a function of increasing temperature and 
can thereby become much smaller than the mean spacing 
between impurities. 



IV. QUANTUM BOLTZMANN EQUATION 

For a fermionic quantum many-body system under 
non-equilibrium conditions, we define the usual Keldysh 
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propagators 



-i{'ip{xi)i>Hx2)) 
i{tlj\xi)ip {X2)) 

e{h-h)G>{xi,x2) 
e{h^t,)G<{x,,x2) 
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+ e{h-t2)G<{xi,X2) 

and formulate a Dyson equation 

G = Go ^" GqYjG 
for the matrix Green's function 
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G 



Gt 
G> 



-G< 
'G-t 



(13) 
(14) 

(15) 



and S defined in analogy. Gq contains the non- 
interacting propagators. We employed the combined 
variable x — (r, t) that parametrizes the fermionic field 
Without further complications, it can also include 
spin and other variables. For brevity, we dropped these 
coordinates in Eqs (14) and (15) as well as we dropped 



integrating over intermediate coordinates in Eq. (14). 
For transport processes, it is convenient to introduce a 
center-of-mass coordinate system 



(R,T) 



1 



{xi + X2) 



{r,t) ^ X1-X2 . 



(16) 



Starting from the Dyson equation ( 14 ) it is then possible 
to derive an equation of motion for G^{xi,X2), which 
we refer to as quantum Boltzmann equation. This has 
been discussed in detail by Mahan.'22l For these purposes, 
it is convenient to work in the Fourier space of relative 
coordinates r 

G(k,w,R,r)= / d'^re'^'' ( dte"^'G{r,t,R,T) . (17) 



A. Derivation 

Our transport description starts from the gauge invari- 
ant quantum Boltzmann equation'^'^ 
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G< - ieE 



-VkS< 
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I]>G< - E<G> 



dijj 



VkReGr, 



+ ^[Re[Eret],G<]+^[S<,Re[G,et]] ■ 
Here, we employed the Poisson bracket 

assG _ (mdG_ 
VkS - VrG + VrS - VkG 



(18) 



P,G] 



and used the retarded propagator Gict defined by the re- 
lation Grot = Gt—G^. Further simplifications to Eq. (18) 
can be made in absence of an external electrical field 
(E = 0), such that all terms proportional to E can be 
dropped. In addition we linearize both hand sides in 
Vr. Here, we just state the result that we justify further 
in Appendix [K\ 

df VrT 



A- 



duj T 



vkcj = E^G"^ 



I]<G> 



(20) 



In the driving terms, the quantities F = — 2ImSi.ct, 
A = i{G> - G<) and the Fermi function /(w) = (1 + 
exp(/3a;))~^, /? = Ijk^T, appear. From the solution to 
this equation, the thermal current density is obtained as 



Q 



(27r)3 TO/ J 27r 



duj 



a;G<(k,a;) 



(21) 



By linearizing this equation in the driving field VrT, it 
is possible to obtain the thermal conductivity k if the 
boundary condition J = is obeyed. 



In order to specify the collision term in Eq. ( 20 1 , we con- 
sider the non-equilibrium self-energy contributions from 
single gauge boson exchange processes 



S^(k,^) 



q ' 



dv 


k X q 


2 






ImZ?(q, v) 


/o I" 


ruf 



{n{y) + l} 



G<(k -I- q, w + z/) + n(j/)G<(k + c\,uj~v) 
and the contribution from potential scattering 



Shnp(k,^^) 



(27r) 



:[rkp(c.)]^G<(p,^) 



(22) 



(23) 



The total self-energy is given by the sum S< = + 
^inip- note that in this approximation, we have ap- 
proximated the gauge field propagator by its equilibrium 
form. This approximation implies that the presence of 
a thermal current of gauge bosons does not change the 
spinon thermal conductivity. No contradiction with pre- 
vious approximations arises. Although without impuri- 
ties, the system has the translational symmetries of the 
Bravais lattice, the finite spinon thermal conductivity we 
will obtain within the approximation of an equilibrated 
gauge field is not artificial. A well-known analogy is 
the electrical conductivity of a Fermi liquid, which has a 
well-established temperature dependence that can be ob- 
tained from scattering due to quasiparticle interactions, 
without involving mechanisms that break translational 
invariance. We will comment further on this approxi- 
mation when we discuss the thermal conductivity of the 
gauge bosons. 



B. Variational principle 



It turns out that Equation (20) is too complicated to 



(19) solve in this form. We attempt to simplify it further such 
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that it can be rewritten in form of a variational princi- 
ple for the distribution function C^Ck^oj). We make the 
linearized ansatz 



/(^) ^ • VkA(k,a;) 



function A(k,ti;), which is sharply peaked as a function 
of ^. This variable can therefore be eliminated from the 
quantum Boltzmann equation, and we can make use of 
the reduced distribution function 



G<(k,w) = iA{k,uj) 

ULU ± 

(24) 

which is valid in the regime |VrT/T| ^ 1. All non- 
equilibrium properties are encoded in the vertex distribu- 
tion function A(k, w) which still needs to be determined. 
First, we assume that the collision term is dominated by 
scattering on the gauge field, ^ ^iiiip- This approxi- 
mation is valid in the temperature regime Ti <C T < T* 
discussed above. Employing Eqs (22) and (24 1, the col- 
lision term can be specified as 



i[Y.<A^ 2TG<] = 2r(k, w)A(k, w 

A(k,^)^ 
q 

df{uj 





k X q 


2 






ImD(q, v) 


la 





df VrT 

du T 

T 



VkA(k,tj) 



Vk4 



l-f{LO + v) 



A(k + q, a; + i/) A(k + q,uj + i^)- [n(-i^) + /{uj-f 
1 - /(^ - i^) 



1 - /(^) 



-y4(k + q,uj - i^)A{k + q, w - ly)] 



(25) 





k X q 


Jo TT 





Solving this equation for the vertex distribution function 
yields the self-consistency condition 

^^^'") = 2f(^;^-2r(^^ 

ImDic, .) ^-^±^ { [.(.) + /(. + .)] ^i^^ 
^(k + q, w + i/)A(k + q, w + I/) - [n{-v) + /(w - v)] 



1 - /(^) 



- A(k -I- q, w - I/) A(k -I- q, w - !^) 



(26) 



Equation |26] is solved by a vertex distribution function 
A(a;) that depends only on frequency but not on mo- 
mentum. This result is of central importance in order to 
formulate a variational ansatz for (k, uj) later on. To 
begin with, r(k, w) does not depend on momentum ac- 
cording to Eqs ([t]) and ([s]), such that A(k-|-q, lo) depends 
on k only through the variable ^ = fk+q — The col- 
lision integral is dominated by gauge boson frequencies 
V ~ q^Xd/^ , with the frequency being set by ~ fcsT. In 
the limit ksT <^ kpXd/^ — ^p, we can therefore approx- 
imate g <C fci?, such that all other terms in the integrand 
of Eq. ( 26 1 approximately depend only on q but not on k. 



with corrections of higher order in the small parameter 
q. We will give an estimate for the spinon Fermi en- 
ergy Cp below. Within this approximation, the integral 
can be entirely parametrized in terms of the variables 
f — ek+q — H and q, and the value of the integral will not 
depend on k. Therefore, G^(k, w) will depend on the 
variable ^ = Ck ^ M only through the equilibrium spectral 



dC 
2tt 



G<(e,k, 



(27) 



In this case, it is possible to represent the thermal con- 
ductivity as a variational functional.'^Sl A similar ap- 
proach has been employed previously by Nave and Lee^ 
in the case of two dimensions. We repeat their steps 
wherever it is instructive in order to understand our re- 
sults for the three dimensional case. 

Introducing an ansatz of the form /(k, w) — f{uj) — 



(k, aj)S^ immediately leads to 



(k,c.) 



VrT 

T 



VkA(a;) 



(28) 



Assuming a power-law dependence A(w) oc w'^j'^the trial 
function is given by 0(k, w) = a;k(VRT/T). Upon in- 
serting the definitions for S and Q into the general form 
of Eq. ( 30 1 , it is seen that the power s in the ansatz 
A(a;) cx drops out of the temperature dependence of 
the thermal conductivity once frequency is rescaled with 
temperature. Therefore, we dropped the exponent s from 
our variational ansatz. Note that it is not necessary to 
normalize '/>(k, uj) since the normalization constant will 
drop out of our final results. We next define the rate of 
entropy density production as 



5" = V 



Q 



(29) 



generated by the thermal current density Q. 
the variational functional 



I mmimizes 



W : 



(30) 



where Q is the heat current density given by 
Q = iV(0) / dkdu}(l)VFUJ^k 

J OUJ 

and the rate S of entropy density production is 

S = —5^-/32 / dujduj'di^dqdkdk'lk' x qP 
mj J 

ImDiq, j.)(0 - 0')'/o(w)[l - foioj')]noiiy) 



(31) 



X 5{uj' - UJ - v)&{{k! -k) - ^) 

kp 



(32) 



Here the integration over k is defined as integration over 
the d — 1 dimensional unit sphere, A''(0) is the den- 
sity of states at the Fermi level for up and down spins 
combined. In 2d, we have N{0) — rnf/n and in 3d 
A^(0) = nifkp / {211^) The global minimum of the vari- 
ational functional W is identical to the spinon thermal 
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resistivity k1^. 

As discussed in Section III C at lowest temperatures the 
energy relaxation rate is impurity dominated 



r(k,u;)^2ro 



(33) 



assuming a constant impurity contribution 2ro — 
2niIm[TkpkF('^ — 0)]- this case, the rate of entropy 
density production in Eq. ( 32 ) has to be modified accord- 



ingly. Instead, this case can be treated in a much simpler 
way by using the Wiedemann-Franz law for free fermions 



(34) 



with the Lorentz number Lq — {n'^ /3){kB /e)'^ ■ We can 
use the Drude formula for the residual conductivity of a 
Fermi gas 



Kof 



(35) 



and insert it into Eq. (34) to determine the spinon heat 
conductivity Kf. The spinon density nf \s equal to the 
inverse unit cell volume, while the spinon scattering rate 
tJ^ can be equated to twice the self-energy correction 
from impurity scattering, 2ro. Independent of dimen- 
sion, up to a number of C(l), the heat current density is 
given by Q = N{0)vF{kBTf- 



In order to determine S in Eq. ( 32 1 , we first integrate over 
the variables lOjUj' and q and parametrize the integrals 
over k and k' in spherical coordinates. Furthermore, we 
can approximate 



(36) 



where v is given by the relative frequency v = uj — uj' . 
Since the main contribution to the integral comes from 
q ^]^g term i^^(k • u)^ dominates the temperature 

dependence at ksT <iC- tp and we can neglect the term 
2v{uj — iif'0{q) and all higher corrections. A straightfor- 
ward but somewhat tedious integration shows then that 
the rate of entropy production follows S ^ (kBT)'^^^'''^^^ . 
We note that we did not introduce an infrared cutoff fre- 
quency for the gauge field mode in Eq. ( 32 1 , although 



this leads to a divergent retarded self-energy at finite 
temperatures.^ Even at one loop order, this divergence 
is canceled in the collision integral of the QBE. Intro- 
ducing an infrared cutoff given by = 1 would only 
renormalize the prefactor of the overall result by a num- 
ber of 0(1), as we checked explicitly. 



C. Thermal conductivity: results 

Putting together the Drude result and the result for 
scattering on gauge bosons, the spinon conductivity is 





(I) T < Ti 


(II) Ti < T < Ta 


(III) Ta < r 


2d 
3d 


ciT 
diT 


d2 


non-universal 
non-universal 



TABLE I: Temperature-dependence of the total thermal con- 
ductivity K in different temperature regimes. Estimates for 
the characteristic temperatures Ti and T2 as well as expres- 
sions for the parameters ci, C2, di and da are discussed in the 
main text. The temperature-dependence in regime (III) is 
non-universal since it depends also on the temperature scale 
T* , among other aspects. It is discussed in more detail in the 
main text. 



given by 



T 





,2/3 3 mt 


^ B 


1 kl 








-1 


h ksT 


3 rri; 




1^1 







(2rf) 
(3d) 



(37) 



Here, we restored units of h and ks and omitted prefac- 
tors of 0(1). The spinon density n/ corresponds sim- 
ply to the inverse of the unit cell volume, such that 
rif = l/Vuc- We defined the unit cell volume Vuc as 
the volume of the two-dimensional unit cell for the two- 
dimensional thermal conductivity and as the volume of 
the three-dimensional unit cell for the three-dimensional 
thermal conductivity. Therefore, we defined the two- 
dimensional thermal conductivity as the thermal conduc- 
tivity of a single two-dimensional layer. It turns out from 
the results in Eq. ( 37 1 that in both two and three dimen- 
sions, the spinon thermal conductivity agrees with the 
basic argument n ^ CyVpTE stated by Nave and Lee,'22l 
assuming that the specific heat follows Cy ~ jT and 
te is the energy relaxation rate of the spinons, which 
is constant for isotropic impurity scattering and follows 
Te ~ y-t^/s for scattering on gauge field fiuctuations.'^ 
The result for the two-dimensional thermal conductivity 
has been already obtained by Nave and Lee.'^S For the 
case of three-dimensional candidate materials, we have 
now obtained the finding that the thermal conductivity 
of the spinons due to scattering on gauge field fluctua- 
tions has a temperature- independent contribution. This 
can be seen from the result for 3d in Eq. (37), since 
Kf/T ~ 1/T in the limit t7^ = 0. However, due to 
impurity scattering the total thermal conductivity will 
vanish in the limit of zero temperature, such that the 
constant term d2 appearing in Table [T] can only appear 
above the temperature scale Ti introduced in[Tj It is im- 
portant to recall that we have approximated the collision 
integral by the sum of impurity and gauge field collision 
processes, and it has been assumed that any interference 
between these types of collision processes does not influ- 
ence the temperature dependence of the spinon thermal 
conductivity. 

Finally, the contribution of conduction electrons to 
the thermal conductivity can be obtained by using the 
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Wiedemann- Franz law ( 34 ) 



(38) 



with the Fermi liquid resistivitjE^ 



(39) 



Here, a is the interlayer spacing of the metallic system, 
ep the conduction electron Fermi energy and the resid- 
ual resistivity po can be determined from the Drude for- 
mula (35). In metallic systems, the Lorentz number Lq 



is up to a factor of 0(1) equivalent to the free electron 
value used in p4|). 



We now summarize the behavior of the total thermal 
conductivity k given by 



K = Kc + Kf ■ 



(40) 



Based on our results, we distinguish the quantitative be- 
havior of three different temperature regimes (I)- (III) 
summarized in Table [H 

Regime (I) is the impurity dominated regime with a lin- 
ear temperature dependence as in a conventional Fermi 
liquid. The ratio of spinon vs. conduction electron con- 
tribution in this regime is 



Kc ruf rifTf 



(41) 



mc/mf is a small factor that is approximately given by 
the ratio Jh/D of the spinon dispersion to conduction 
electron band width D. Here, Jh is understood as the 
strength of nearest neighbor RKKY exchange that is nat- 
urally small compared against the conduction electron 
band width. This could yield to a dominance of the con- 
duction electron contribution to the impurity dominated 
thermal conductivity, although the relaxation times also 
play a role in this ratio. 

Regime (II) is the non-Fermi liquid regime dominated 
by the temperature dependence of the spinon contribu- 
tion which extends from Ti to some larger temperature 
scale T2- It shows a rather sluggish T'^l'^ dependence in 
two dimensions and a constant contribution in three di- 
mensions. The temperature-dependence in regime (HI) 
depends on various degrees of freedom, such that the 
temperature scale T2 depends not only on parameters of 
the spinon system. First, the leading temperature depen- 
dence of the spinon thermal conductivity in regime (II) 
is only valid at temperatures fc^T < Cp^ as we explained 
in the context of equation (26). Next, spinon scattering 



on hybridization bosons can only be neglected below the 
temperature scale T* . We note that also gauge bosons 
carry a finite thermal current that leads to corrections to 
the total thermal conductivity. As discussed below, this 
correction is again negligible for temperatures fc^T < e^p. 



We can estimate the spinon Fermi energy tp by equat- 
ing it to the typical energy scale Jh for RKKY exchange 
interaction, tp w Jh ~ lOOiiT. Finally, the temperature 
dependence of the conduction electron contribution has 
to be compared against the temperature dependence of 
the spinon thermal conductivity. We use the estimate 



(42) 



knT 



for the temperature scale when both contributions be- 
come equal and obtain fc^T = £'p{4/£p)'''^ Since 
(ip/ip) <C 1, this temperature scale is much larger than 
the conduction electron Fermi energy, which is of order 
10*K, and we can therefore neglect the temperature de- 
pendence of the conduction electron thermal conductivity 
in this case. We conclude that the temperature scale T2 
is given by the minimum of the energy scales ep and T* . 
A more detailed discussion of regime (HI) relies on addi- 
tional high energy degrees of freedom that we excluded 
from our theory so far, and a more detailed analysis of 
this temperature regime is required to accurately discuss 
the thermal conductivity beyond the temperature scale 
T2. 



D. Wiedemann- Pranz ratio 

A suitable quantity to detect spinon contributions to 
the thermal conductivity is the Wiedemann-Franz ratio 



L{T) = 



lit 

<JtT 



(43) 



Due to the absence of a spinon charge current, the 
presence of a spinon thermal current will always imply 
L{T) > Kc/{(JcT). Since in our theory the conduction 
electron Wiedemann-Franz ratio Kc/{o'cT) behaves as in 
a Fermi liquid, the additional spinon contribution will 
lead to characteristic deviations from Fermi liquid be- 
havior. In the following, we distinguish two characteristic 
regimes for the Wiedemann-Franz ratio that are identical 
to the regimes (I) and (II) for the thermal conductivity 
of the last subsection. 

In regime (I), the impurity dominated regime, the 
Wiedemann-Franz ratio is given by 



<JtT 



Ml + ^) 



(44) 



with Uf — {nfTf)/mf and Uc — {ncTc)/mc. The absence 
of a spinon electrical current causes a renormalization of 
the Wiedemann-Franz ratio in the impurity dominated 
regime that might be unobservably small, due to the 
small mass ratio rric/mf ^ Jh/D <C 1. 
More promising with respect to observability is regime 
(II), where spinon collisions are dominated by scatter- 
ing on gauge bosons. There, the Wiedemann-Franz ra- 
tio is given by the ratio of temperature-dependent parts 
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FIG. 2: Sketch of the characteristic temperature depen- 
dence of the Wiedemann- Franz ratio L{T) in the temperature 
regimes (I) and (II) discussed in the main text. The quantity 
L{T) corresponds to the full line. At zero temperature, the 
ratio saturates at the finite value L{T — 0) = Lo{l + -^). 
At a temperature Tmax that is i.e. dependent on the spinon 
and conduction electron scattering rate, a maximum occurs 
that turns over into a characteristic power law divergence 
y-(6+d)/3 ^y^^^ sketched in form of the dashed hue. 



of spinon thermal conductivity and conduction electron 
electrical conductivity. In this case the Wiedemann- 
Franz ratio shows the temperature dependence 



LiT) 




(45) 



These divergences will be cut off at low temperatures 
by the thermal conductivity due to impurity scattering. 
Therefore, the behavior stated in Eq. (45 1 leads to a char- 



acteristic maximum of the Wiedemann-Franz ratio at a 
temperature T^ax that can be estimated by equating the 
two contributions to Eq. (37). This estimate results in 
the temperature 



B ^ max 



HUfTf 



in two spatial dimensions and 



hnifTf 

in three spatial dimensions. More general, any thermal 
conductivity with a temperature dependence dominating 
the Fermi liquid T'^-dependence will lead to a maximum 
of the Wiedemann-Franz ratio. This assumes that the 
electrical conductivity of the system behaves like in a 
Fermi liquid. 



E. Corrections to thermal conductivity 

The contribution of gauge bosons to thermal transport 
can be analyzed in analogy to the spinon transport equa- 



tion by defining a distribution function 

£D<{q,iy,r,t) . (46) 

It is now again possible to use the variational princi- 
ple used already in this section to calculate the spinon 
thermal conductivity if the spinous are assumed to be in 
equilibrium in the gauge boson transport equation. How- 
ever, in principle the gauge boson thermal current leads 
to a modification of the spinon thermal current and vice 
versa. This effect could be included by considering the 
coupled system of transport equations for spinous and 
gauge bosons. In analogy to the phonon drag in met- 
als, such corrections to the spinon thermal conductivity 
can be shown to be of relative size {kBT)/ep and can be 
neglected within this accuracy.l^^lSSI jg j^qw completely 
analogous to the calculation of the spinon thermal con- 
ductivity to rewrite the gauge boson thermal conductiv- 
ity in form of a variational functional. We do not repeat 
these steps in detail here, but cite just the result we ob- 
tained. After solving the integrals for the gauge boson 
heat current density and entropy production, the gauge 
boson thermal conductivity Kg can be estimated as 



kuT 



-1 (d+l)/3 



(47) 



with a non-universal constant Cg that depends on details 
of the spinon dispersion. For temperatures ^ 1, the 

gauge boson thermal conductivity Kg is therefore only a 
subleading correction to the thermal conductivity of the 



spinous, see Eq. (37). 



Finally, phonon contributions to the thermal conductiv- 
ity are proportional to [T/Q^Y and are therefore negli- 
gible for temperatures much smaller than the Debye tem- 
perature Qd- 



V. MAGNETIC SUSCEPTIBILITY 
A. Results 

In our model, the physical spin susceptibility is given 
by the sum of the spinon susceptibility and the conduc- 
tion electron susceptibility inside the FL* phase. 



Xs = Xc + Xf 



(48) 



In analogy to our discussion of the thermal conductivity, 
corrections to this formula arise due to the presence of hy- 
bridization bosons and are therefore again exponentially 
small as a function of temperature below the temperature 
scale T*. We will neglect these effects in the following, 
considering again only temperatures T < T* . Part of the 
behavior of the spin susceptibility Xs has been discussed 
already by Paul et aLj^iwith a focus mainly on the finite 
temperature regime above the QCP. In the FL* regime 
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considered by us, the conduction electron contribution 
Xc is described by the usual Fermi liquid form 



2 2A^c(0) 



The paramagnetic part is controlled by the conduction 
electron density of states Nc{0) at the Fermi level and 
the Landau parameter Fg. Different from a Fermi gas, 
in two dimensions linear temperature corrections arise in 
a Fermi liquid,'^ while in three dimensions, the usual T^ 
dependence prevailsP^ The coefficients bi and 62 depend 
on interactions and band structure details and while hi 
is positive,'^ the sign of 62 has not yet been found to 
the best of our knowledge. One way to obtain the spinon 
contribution Xf is the free energy of the system of spinons 
and gauge field 



F ^ -ksTlnZ Z = Va Vfe 



(49) 



In RPA approximation, this free energy is given by + 
Fa, where Fa is the correction to the free fermion part 
Fs due to quadratic fluctuations of the transverse gauge 
field, 



(2^ 



In 



(50) 



Therefore, Xf "wiH have the form Xf = X/ + ^X/j where 
X° — {gj /2)^\N f{0) has the Fermi gas form of free 
spinons with g- factor g f and density of states at the Fermi 
level Nf{0). The contribution Axf denotes additional 
fluctuation corrections to the susceptibility. 

It yet remains to discuss the dependence of Fa on 
an external magnetic field. In presence of an external 
magnetic field, Fa will split into two contributions Faa 
stemming from different spin projections a = zt^. The 
magnetic ffeld couples to the spinons by a Zeeman term 
Fz = gffJ-nB that splits the Fermi wavevector to the 
spin-dependent values kpc ~ \/kjp + a2mfEz, depend- 
ing on the spinon g- factor gj. In this way, it is straight- 
forward to obtain the shift Axf from 



Ax/ 



d^Fa 



(51) 



3=0 



The low temperature behavior is 



Ax/ 




aiT'i (2d) 
a2T2log(r/e;^) (3d) 



Details of the calculation are shown in App. [B] where in 
particular the form of the parameters ai and 02 is clar- 
ified. Both fli and 02 are positive and depend on the 
spinon dispersion. The temperature-independent term 
Xo is a- paramagnetic contribution that has been dis- 
cussed in detail in Ref. [231 It is roughly of the same 
order of magnitude as the Pauli paramagnetic suscepti- 
bihty x°f 



B. Discussion 

Possibly the best way to observe the spinon contribu- 
tion to the susceptibility is to consider the zero temper- 
ature limit. In this limit, the ratio of the spinon to con- 
duction electron Pauli susceptibility is given by the ratio 
of their densities of states at the Fermi surface. This ra- 
tio can be estimated from the ratio of their band widths, 
D/Jh- The spinon band width energy scale Jr is gen- 
erated by RKKY exchange according to Jh ~ ^kZ-D, 
such that the spinon Pauli susceptibility dominates that 
of the conduction electrons by a large factor / ^ 1. 
This strong effect could therefore be used for the exper- 
imental identification of such an FL* phase with gapless 
fermionic spinons. Indeed, in the non-Fermi liquid phases 
recently observed in experiments on Yb-based materi- 
als, the spin susceptibility at lowest applied temperature 
was enhanced by up to a factor of thre^ and even fiv^ 
as compared to the Pauli susceptibility deep inside the 
heavy Fermi liquid phase. 

The temperature-dependent part of the spinon suscep- 
tibility also leads to important modifications of the full 
temperature dependence of the magnetic susceptibility. 
In two dimensions, the asymptotic temperature depen- 
dence in the limit T — )■ is set by the conduction electron 
contribution and behaves linear in temperature. The free 
spinon part X/ does not to contribute to this asymptotic 
behavior, since it has the characteristic quadratic tem- 
perature dependence of band fermions. Still, it will in- 
fluence the overall temperature dependence of the spin 
susceptibility. Since in general the spinon Fermi en- 
ergy is quite small compared to the conduction elec- 
tron Fermi energy, the spinon contribution Ax/ might 
qualitatively change the temperature dependence of the 
susceptibility at quite low temperatures. An estimate 
is that the Ts -dependence due to gauge field correc- 
tions is dominating above a temperature scale set by 
ksT/e'p = [kBT/e^pf/^, such that 



J 



f \ 3/2 



(52) 



Estimating tp ~ Jh ~ lOQK and w IQ^K we get 
T sa Q.IK for the crossover to the spinon dominated 
temperature-dependence of the spin susceptibility. 
For three-dimensional candidate materials, both the tem- 
perature dependence of the gauge field contribution and 
the band susceptibility of the spinons should dominate 
the temperature dependence of the conduction electron 
contribution to the total spin susceptibility. The contri- 
bution Xj is expected to do so because it is proportional 

{ksT /ep)"^ , with the lower Fermi energy ep <^ Cp. Ul- 
timately, the logarithmic enhancement ~ \og{kBT/ep) 
set by the gauge field corrections is expected to dom- 
inate the temperature dependence of Xs in the regime 
kgT/ep ^1. In general, our results predict a positive 
finite temperature correction to the Pauli susceptibility. 
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VI. 



SUMMARY AND DISCUSSION 



ground state is destroyed. 



We have discussed the thermal conductivity and spin 
susceptibility of a U(l) fractionalized Fermi liquid in 
two and three spatial dimensions. Within the scenario 
of gapless fermionic spinon excitations, the total ther- 
mal conductivity becomes enhanced by the spinon con- 
tribution already in the impurity dominated regime, al- 
though the enhancement might be unobservably small 
due to the small Fermi velocity of the spinous. In a 
temperature regime where scattering of spinous is domi- 
nated by gauge boson absorption and emission, we were 
able to argue that the spinon thermal current will domi- 
nate the thermal conductivity as compared to conduc- 
tion electron, phonon and gauge boson contributions. 
For the important three-dimensional candidate state, 
we obtained a characteristic plateau in the tempera- 
ture dependence of the heat conductivity that dominates 
all other temperature-dependent contributions in mag- 
nitude. We propose to detect such anomalous contri- 
butions to the thermal conductivity by measuring the 
Wiedcmann-Franz ratio. As a function of temperature, 
this ratio shows a pronounced temperature maximum 
which is caused by the non-Fermi liquid thermal con- 
ductivity of the spinon system. The size of the max- 
imum depends strongly on the impurity concentration 
and can largely exceed the Wiedemann-Franz ratio of 
conventional metals if the system is sufficiently clean. 

Our theoretical predictions particularly rely on the ex- 
istence of gapless fermionic excitations. However, many 
other types of spin liquid states might be realized in na- 
ture with gapped/and/or bosonic excitation spectra that 
do not show a Pauli paramagnetic susceptibility. From 
our analysis of the spinon susceptibility in the U(l) FL* 
phase, we predict a large enhancement of the Pauli para- 
magnetic susceptibility due to the flat dispersion of the 
fermionic spinon excitations caused by RKKY exchange. 
This enhancement should be observable at temperatures 
ksT < ep. Such a large paramagnetic contribution is un- 
likely to be observed in presence of bosonic excitations 
or a gapped spectrum. A gapped spinon spectrum would 
be characteristic for an increasing spin susceptibility as 
a function of temperature if the size of the gap is smaller 
than the spinon Fermi energy.'22! Such a behavior is not 
observed in recent experiments.'^ 

As an additional thermodynamic signature for our 
three dimensional candidate state, however, we predict a 
logarithmic enhancement of the temperature-dependent 
part of the Fermi liquid susceptibility from the spinon 
contribution. Our results might be complemented by 
future work that may help to predict other signatures 
of ordered states in the local moment system coexisting 
with itinerant conduction electrons. A unique probe for 
fermionic spin liquids may be the thermal Hall effect 
Undoubtedly, this demands further experimental effort. 
Moreover, microscopic calculations need to be established 
that can explain the suppression of broken symmetries 
in Kondo lattice systems once the heavy Fermi liquid 
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Appendix A: Derivation of quantum Boltzmann 
equation 



It is possible to rewrite the collision term in ( 20 ) as 
I]>G< - I]<G> = -i[2TG< - S<A] . 



(Al) 



According to our assumption, the gauge field propagator 
is replaced by its equilibrium form, and the drag between 
spinous and gauge field is neglected. Furthermore, the re- 
tarded spinon propagator obeys the equation of motiorP^l 



UJ ■ 



8m/ 



V R 



Gr, 



1 . 



(A2) 



In the linearized spatially homogeneous quantum Boltz- 
mann equation, the retarded self-energy Ei-ct(k,w) can 
therefore be replaced by its equilibrium form, since field 
strength enters only to ©(E^) or 0(V|{^). Since the col- 
lision term (Al) has to vanish in equilibrium, only field 



dependent terms have to be retained, and the collision 
term reads therefore 



i[2TG< - ^<A] = 2T{k,uj)A{k,uj) 



du! 



F-VkA(k,w) 



+ A{k,Lu) 



E 

■ g 



dv 
vr 



k X q 



nif 



ImZ3(q, //) 



Vk+q 



1 - /(^ + 



1 - /(^) 

A{k + q,uj + v)K{k + q,Lj + v) - [n{-v) + /(a; - v)] 
1 - /(^) 



-A{k + q,u}- v)K{k + q,uj -v)] \ . (A3) 



The driving force F is given by the ratio VT /T in our 
concrete case. It remains to derive the driving term of 
the QBE (pOl). 



In presence of a thermal gradient, we can still neglect any 
dependence of the QBE on the absolute time variable T. 
However, the spatially dependent temperature distribu- 
tion causes a dependence of the functions and G^ on 



13 



the coordinate R. The driving term therefore includes 
the terms 

- VRRcErct • VkG< + Vk(ek + ReS.ct) • VrG< 

- [S<,ReGrct] . (A4) 

We can use the arguments of Ref. [2H in order to show 
that the term [S^, ReGjct] has no influence on physical 
observables. Moreover, the gradients VRReErct and 
VkReSrot vanish according to the assumptions stated 
above. Therefore, the driving term is simply given by 
VkekVrG<. 



In three dimensions, we get 
FiT) F(0) = ^ (^) ) + OiT') (B3) 

with a momentum cutoff A that is of order the spinon 
Fermi wavevector kp. The quantity XdA'^/F is equivalent 
to the spinon Fermi energy up to a factor of 2/tt. Here 
and in the following, we employ units where fcs = 1. In 
a finite magnetic field, the ratio depends on the spin 
projection a, 



Appendix B: Calculation of susceptibility 

Following the same steps as in Ref. IH we obtain for 
the temperature-dependent part of the free energy 



4: f d'^k r 
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The effect of a magnetic field can be analyzed afterwards 
by considering the field-dependence of the parameters F 
and Xd- 

We introduce the new variables x — and u = ^ and 
get in 2d 



F{T) - F(0) = Ta 



F y F(|,0) 
Xd J 37r2 



13.3425 -hC'(r^p2) 



£o_ _ I Unkpamf {2d) 
Xda \ nkp^nif {3d) 



with the shifted 



Fermi wavevectors kp^r = 
The temperature-dependent 
part of the susceptibility (51) is finally obtained by 
using 



agf2mf^BB. 
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The temperature-independent contribution to Eq. (51) 
can be obtained in a similar way^^- 
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